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1 INTRODUCTION 



ABSTRACT 

We study the cosmological information of weak lensing (WL) peaks, focusing on two 
other statistics besides their abundance: the stacked tangential-shear profiles and the 
peak-peak correlation function. We use a large ensemble of simulated WL maps with 
survey specifications relevant to future missions like Euclid and LSST, to measure and 
examine the three peak probes. We find that the auto-correlation function of peaks 
with high signal-to-noise (S/Af) ratio measured from fields of size 144 deg 2 has a 
maximum of ~ 0.3 at an angular scale i9 ~ 10 arcmin. For peaks with smaller S/Af, 
the amplitude of the correlation function decreases, and its maximum occurs on smaller 
angular scales. The stacked tangential-shear profiles of the peaks also increase with 
their S/Af. We compare the peak observables measured with and without shape noise 
and find that for S/Af ~ 3 only ~ 5% of the peaks are due to large-scale structures, 
the rest being generated by shape noise. The correlation function of these small peaks 
is therefore very weak compared to that of small peaks measured from noise-free 
maps, and also their mean tangential-shear profile is a factor of a few smaller than 
the noise-free one. The covariance matrix of the probes is examined: the correlation 
function is only weakly covariant on scales i? < 30 arcmin, and slightly more on 
larger scales; the shear profiles are very correlated for i? > 2 arcmin, with a correlation 
coefficient as high as 0.7. The cross-covariance of the three probes is relatively weak: the 
peak abundance and profiles have the largest correlation coefficient ~ 0.3. Using the 
Fisher-matrix formalism, we compute the cosmological constraints for {f2 m , erg, u>, n s } 
considering each probe separately, as well as in combination. We find that the peak- 
peak correlation and shear profiles yield marginalized errors which are larger by a 
factor of 2 — 4 for {O m , erg} than the errors yielded by the peak abundance alone, while 
the errors for {w, n s } are similar. By combining the three probes, the marginalized 
constraints are tightened by a factor of ~ 2 compared to the peak abundance alone, 
the least contributor to the error reduction being the correlation function. This work 
therefore recommends that future WL surveys use shear peaks beyond their abundance 
in order to constrain the cosmological model. 



weak gravitational 
Euclid and LSST 



With the approach of large 
lensing (WL) missions, such a s 
JLSST Science Collaborations et all 120091 ; lLaureiis et al.l 
l201ll ). it is imperative to optimize the extraction of cosmo- 
logical information from WL probes. Until now, the WL 
community has been using mainly the 2-point function of 
the shear field in order to obtain cosmological constraints 
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(Jarvi s et al.l 12 003: Hockstra et al.l 120 06: Scmboloni et al.l 
120061 : iHetterscheidt et alj|2007l : iKilbinger et al.ll2012fl. and 
also to plan and forecast coming missions l|Laureiis et al.l 
l201ll h Whilst a lot of effort is concentrated on scrutinizing 
the systematics of WL surveys - a most essential endeavour 
for the success of the next-generation missions, and indeed 
for the future of WL as a cosmology probe - the task of se- 



lecting the most efficient statistics for parameter estimation 
is also a crucial aspect that should not be neglected. 

One of the possible WL observables involves the 'peaks' 
of the shear field, i.e. regions of high signal-to-noise (S/Af) 
of the field, produced by overdense regions of the den- 
sity field projected along the line of sight. Shear peaks are 
the imprint of clusters in WL maps, and can be used to 
detect and measure cluster masses llHamana et al. 2004 
" 20041: iMaturi et all [20051: iHeimawi fc Spergel 
Marian fc Bernstein 2006; Schirmer et al 



Wang et al 

20051: Dahlell2006l _ 

20071 ; Maturi et al.ll2007l : lAbate et al.ll2009f h The peak abun 



dance scales with cosmological pa rameters in the same 
way as the halo mass function (|Reblinskv et al . 1999; 
M arian et all 120091 . l2010h . and therefore can be used to 
constrain the cosmological model l|Dietrich fc Hartlapll2010l ; 
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iKratochvil et alj|2010l ; iMarian et ai1l2012l ; iBard et alj2013j ) . 
The shear-peak a bundance can also c onstrain primordial 
non-Gaussia nitv (IMarian et alJ l201ll ; iMaturi et all l201ll : 
Hilbert et al ]|2012l ), and so far seems to be possibly the most 
effective WL observable suitable for that purpose - although 
more work must be done on this subject. 

Whilst the abundance of peaks has been investigated 
by a number of studies, higher-order statistics of the shear 
peaks have until now been overlooked. Yet, since the abun- 
dance of peaks can be used as a cosmological tool, one would 
expect their clustering to also be valuable. In analogy with 
3D observables, just as both the halo correlation function 
and the halo mass function are very sensitive to cosmology, 
the same could be true about the shear peaks. 

The main goal of this study is to investigate the correla- 
tion function of WL peaks, and in particular to quantify the 
improvements that it yields on cosmological constraints rel- 
ative to the peak abundance. As a secondary line of inquiry, 
we shall also pursue the cosmological information contained 
in the tangential-shear profiles of the peaks. While these 
statistics - correlations and stacked profiles of overdensities 
- have been already studied in various the oretical works, 
as well as by using real and simulated data |Sheldon et al.l 
2004 1 2009: Mandel baum et al.|[201(il: lOguri fc Takadall201lF 
Mandelbaum et al]|2012l : ICacciato et al.ll2012T ). here we ex- 
plore the characteristics of shear-selected overdensities for 
which the only known information is the redshift distribu- 
tion of the source galaxies. This might seem a somewhat 
unlikely and conservative scenario for future surveys, but 
here we just examine the intrinsic information contained in 
the peaks, even if our results will be sub-optimal on this 
account. We defer to a near- future study a more realistic 
scenario, where prior information on the peaks is available, 
and where we can account for various survey uncertainties, 
such as intrinsic alignments, masked regions, calibration er- 
rors of the source galaxies, etc. 

We shall present measurements of the peak function, 
peak-peak correlation function, and peak profiles from an 
ensemble of WL maps generated with ray-tracing through 
iV-body simulations with varying cosmologies, as described 
in fj2] We shall use the Fisher-matrix formalism to establish 
the cosmological sensitivity of the above-mentioned probes. 
We underscore that our results are drawn from collaborative 
efforts, which involved generating the iV-body simulations, 
the ray-tracing maps, developing the hierarchical algorithm 
for WL peaks; the present study is part of a larger research 
program, wh i ch has also yiel ded the result s pub lished in 
IMarian et all (|201ll . 120121 ) and lHilbert et all (|2012h . There- 
fore, we shall not give a full description of all our tools, but 
simply refer to these publications. 

The paper is partitioned as follows: in Sj2] we provide a 
description of the ray-tracing maps used in this study; in Ej3] 
we illustrate our measurements of the three peak probes; in 
3D we present the cosmological constraints extracted from 
the three probes; in ^5]we summarize and conclude. 



2 SIMULATED WL MAPS 

We generated WL maps from ray-tracing through iV-body 
simulations. We used 8 simulations which are part of a 
larger suite performed on the zBOX-2 and zBOX-3 super- 



computers at the University of Zurich. For all realizations 
11 snapshots were output between redshifts z = [0,2]; fur- 
ther snapshots are at redshifts z = {3,4,5}. We shall refer 
to these simulations as the zHQRIZQN si mulations, and the y 
were described in detail in lSmithl l|2009h : ISm"ith et all l|2012l ). 
Each of the zHORIZON simulations was performed using the 
publicly available Gadget-2 code (|Springell [20051 ), and fol- 
lowed the nonlinear evolution under gravity of N = 750 3 
equal-mass particles in a comoving cube of length L s i m = 
1500 h~ Mpc; the softening length was i so f t = 60 /i _1 kpc. 
The cosmological model was similar to that d etermined by 
the WMAP experiment (jKomatsu et al.ll2009h . We refer to 
this cosmology as the fiducial model. The transfer func- 
tion for the simulations was generated using the p ublicly 
available cmbfast code (Seli ak fc Zaldarriaeal Il9 96). with 
high sampling of the spatial frequencies on large scales. Ini- 
tial conditions were set at redshift z = 50 using the se- 
rial version of the publicl y available 2LPT code (|Scoccimarrol 
1 19981 : ICrocce et all l2006h . Table [1] summarizes the cosmo- 
logical parameters that we simulated. We also use another 
series of simulations, identical in every way to the fiducial 
model, except that we have varied one of the cosmologi- 
cal parameters by a small amount. For each new set we 
have generated 4 simulations, matching the random real- 
ization of the initial Gaussian field with the correspond- 
ing one from the fiducial model. The four parameter vari- 
ations were: {n -» (0.95, 1.05), ct 8 -» (0.7,0.9), fi m -» 
(0.2,0.3), w — > (—1.2,-0.8)}, and we refer to each of the 
sets as zHORIZON-Vla.b,. . . ,zH0RIZ0N-V4a,b, respectively. 

For the WL si mulations, we consi dered a survey 
similar to Euclid (lLaureiis et alJ l201ll ) and to lsst 
(LSS T~Science Collaborations et al.l 120091 ) . with: an rms 
a 7 = 0.3 for the intrinsic image ellipticity, a source num- 
ber density n = 40arcmin~ 2 , and a redshift distribution of 
source galaxies given by V(z) = N(zo,f$) z 2 exp[— (z/zo) 13 ], 
where the normalization constant J\f ensures that the in- 
tegral of the source distribution over the source redshift is 
unity. If this interval extended to infinity, then the normal- 
ization could be written analytically as: N = 3/(zoF[(3 + 
/?)//?]). There is a small difference between this value and 
what we actually used, due to the fact that we considered 
a source interval of [0,3]. We took f3 = 1.5, and required 
that the median redshift of this distribution be z mc d = 0.9, 
which fixed 20 * 0.64, and gave a mean of z m oan = 0.95. 

From each iV-body simulation we generated 16 indepen- 
dent fields of view. Each field had an area of 12 x 12deg 2 
and was tiled by 4096 2 pixels, yielding an angular resolu- 
tion 6>pi x = 10 arcsec. For each variational model, the to- 
tal area was of 9000 deg 2 , while for the fiducial model it 
was of * 18000 deg 2 . The effective convergence k in each 
pixel was calculated by tracing a light ray back through the 
simulation with a multipl e-lens-plane ray-tracing algorithm 
l|Hilbert et al.ll2007l . [2009) . Gaussian shape noise with vari- 
ance o- 2 /(n# 2 ix ) was then added to each pixel, creating a 
realistic noi se level and correla tion in the filtered conver- 
gence field l|Hilbert et al.ll2007f ). We keep the shape noise 
configuration fixed for each field in different cosmologies, in 
order to minimize its impact on the comparisons of the peak 
abundances measured for each cosmology. 
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Table 1. zHORIZON cosmological parameters. Columns are: density parameters for matter, dark energy and baryons; the equation of state 
parameter for the dark energy; normalization and primordial spectral index of the power spectrum; dimensionless Hubble parameter. 



Cosmological parameters 


a) m 


^DE 


U b 


w 


&8 


n 


Ho[kms _1 Mpc _1 ] 


zHORIZON-I 


0.25 


0.75 


0.04 


-1 


0.8 


1.0 


70.0 


zHORIZON-Vla/Vlb 


0.25 


0.75 


0.04 


-1 


0.8 


0.95/1.05 


70.0 


zH0RIZ0N-V2a/V2b 


0.25 


0.75 


0.04 


-1 


0.7/0.9 


1.0 


70.0 


zH0RIZ0N-V3a/V3b 


0.2/0.3 


0.8/0.7 


0.04 


-1 


0.8 


1.0 


70.0 


zH0RIZ0N-V4a/V4b 


0.25 


0.75 


0.04 


-1.2/-0.8 


0.8 


1.0 


70.0 



3 WL PEAKS AS A COSMOLOGICAL PROBE 

As stated in ijl] we shall consider the cosmological informa- 
tion contained in three shear peak probes: the abundance, 
the profiles, and the correlation function. Until now, only the 
abundance of peaks has been examined as a probe for cos- 
mology; in this study we shall investigate the cosmological 
information contained in the clustering of peaks and their 
profiles, as well as various combinations of the three probes. 

The peaks are dete cted with an ap erture-mass filter, 
i.e. a compensated filter (Schneider 1996), as points of max- 
ima in the smoothed convergence field. They are assigned 
a unique S/M and mass using the hierarchical algorithm. 
This method, as well as the NFW-shaped fi lter function, 
were described in detail in our previous work (|Marian et al.l 
2012). The hierarchical algorithm uses a hierarchy of filters 
of different size, from the largest down to the smallest, to 
determine the size of the NFW filter that best matches each 
peak. For every filter used, we write the aperture mass and 
the S/M at a given point do as 




Sd 2 g[K m (g)-K m (6> A M6> (1 -6>) 

W $ d 2 9 [/c m (0) - « m (0 A )]/t(0o - 0) 



(1) 

(2) 



where M m is the mass of the NFW halo used as a model for 
the filter, and 6a is the aperture radius of the filter, which in 
our case was chosen to be the angular size of its NFW radius. 
K m is the convergence profile of the model, K m is the mean 
convergence inside a certain radius, and K is the measured 
convergence field. If Oo is the location of a peak generated by 
an NFW halo matching the filter, then the above equations 
become 



Af ap (0o) = M a 



S/M(0 o ) = 



(3) 



(4) 



The hierarchical method enhances the detection of peaks, 
their abundance as a function of mass or S/M is independent 
of any particular choice of filter size, and the cosmological 
constraints arising from the peak counts are tighter com- 
pared to the case when a single-sized filter is used. It is also 
instructive to present measurements of the peak function, 
the profiles and the peak-peak correlation function from 
noise-free maps. In that case, the peaks are selected based 
on their aperture-mass values, and they are as signed a hier- 
archic al mass through Eq. as described in lMarian et al.l 
(2012). For the purpose of comparison with the results from 



the noisy maps, we can assign the noise-free peaks a ficti- 
tious S/M value, using EqsQ and @ with the same level 
of shape noise as in the noisy maps. 



3.1 The peak-peak correlation function 

We have measured the correlation function of the peaks 
using the minimum-variance and unbiased estimator intro- 
duced bv lLandv fc Szalavl (|1993f ). The estimator is: 



{DD - 2DR + RR)/RR, 



(5) 



where DD, DR, and RR represent the peak-peak, peak- 
random peak, and random peak-random peak pair counts, 
respectively. The correlation function is measured in 10 an- 
gular bins, logarithmically spaced between t? e [4, 90] ar- 
cmin. The lower bound was chosen to be the smallest possi- 
ble without discreteness effects contaminating the measure- 
ments - the latter appear when the size of the lower bins is 
only a small multiple of the size of the pixels of the map. 
The random catalogue that we created in order to build the 
correlation function estimator contained a million random 



peaks placed in a field of 12 x 12 deg 2 



the same size as 



the simulated maps. 

In Figure[T]we present measurements of the correlation 
function of peaks detected from the maps corresponding to 
the fiducial cosmology. All symbols denote the average of 128 
realizations, the errors being on the mean. The estimator 
of the average correlation function is defined by Eq. (|12p ; 
a similar estimator will also be used to obtain the average 
tangential-shear profile of the peaks and the average peak 
function presented later in this section. 

The left panel of Figure[T] shows the noise- free results, 
while in the right one shape noise is included. The solid 
orange squares/green triangles in the left panel depict the 
auto-correlation cjp p of the peaks with 3.4 ^ S/M ^ 2.6 and 
S/M ^ 3.4, respectively. The solid red circles, blue squares, 
and green triangles in the right panel correspond to peaks 
with S/M > 4.75, 4.75 > S/M > 3.4, and 3.4 > S/M > 
2.6 respectively. The total number of peaks in these bins is 
~ 10, 000, 72, 000, 324, 000 respectively. In both panels, the 
violet crosses depict the cross-correlation Wp p of the bins 
3.4 ^ S/M > 2.6 and S/M > 4.75. 

The correlation functions are strongly negative for the 
smallest bins due to an exclusion effect arising from our peak 
selection: we discard those peaks whose centres are separated 
from the centre of a larger peak by a distance smaller than 
approximately one virial radius. The reason for this choice 
is to avoid classifying substructures as independent peaks. 

In general, the peaks with higher/lower S/M have 
higher/lower correlations, similar to the clustering prop- 
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5 10 50 100 5 10 50 100 

i9]arcmin | i3[arcmin| 



Figure 1. The correlation function of the peaks detected in the maps corresponding to the fiducial cosmological model. The results 
show the average of 128 fields of 12 x 12deg 2 and the error bars represent errors on the mean. Left panel: The orange squares/green 
triangles denote measurements in the absence of shape noise, for peaks with S/Af > 3.4, and 3.4 > S/Af > 2.6 respectively. Right 
panel: Measurements from the noisy maps: the solid red circles, blue squares, green triangles represent the auto-correlation of peaks with 
S/Af > 4.75, 4.75 > S/Af > 3.4, 3.4 > S/Af > 2.6 respectively. In both panels, the purple crosses depict the cross-correlation of the peak 
populations with 3.4 > S/Af > 2.6 and S/Af 5= 4.75, where we have used Eqs ||3} and Q to assign an S/Af value to the peaks. 




0.5 1 5 10 50 0.5 1 5 10 50 

i)[arcmm] i?[arcmin] 



Figure 2. Tangential-shear profiles for the peaks detected in the maps corresponding to the fiducial cosmological model. Left panel: The 
effect of shape noise on the profiles. The red solid pentagons represent peaks detected in the noise-free maps, with S/Af > 4.75. The 
green empty squares depict the profiles of these very same peaks, detected in the noisy maps. Similarly, the orange solid squares show 
the profiles of all the peaks that we detected in the noise-free maps, while the purple stars show the profiles of the same peaks when 
detected and measured in the noisy maps. Right panel: The profiles measured from the noisy maps: the solid red circles, blue squares, 
and green triangles represent peaks with S/Af indicated in the legend. The purple circles depict all the measured peaks binned together. 



erties displayed by 3D halos. For the noiseless measure- 
ments, the maximum of the auto-correlation of peaks with 
S/Af ^ 3.4 is reached at ~ 8 arcmin and is about 0.35, while 
for the peaks with 3.4 ^ S/Af ^ 2.6 the maximum is about 
0.25 and it is shifted towards smaller scales ~ 5 arcmin. The 
auto-correlation functions decrease monotonically to with 



increasing angular scales. The same pattern is followed by 
the correlation functions of the peaks measured from the 
noisy maps: the higher the S/Af bin, the larger the ampli- 
tude of the function and the angular position of its maxi- 
mum. For the largest peaks, a maximum of ~ 0.3 is reached 
at ~ 10 arcmin. Comparing the measurements from noise- 
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free and noisy maps allows us to understand how shape noise 
impacts the correlation function of the peaks. For example, 
the green triangles in the left and right panels of Figure[T]are 
dissimilar, although peaks from the same S/M bin are used 
for the measurements. The noise-free correlation function 
is significantly higher, dropping to only at ~ 100 arcmin, 
unlike its noisy counterpart, which vanishes at ~ 8 arcmin. 
The vanishing owes to the fact that most peaks in this S/M 
bin are generated by shape noise - see also Figure[3]- and 
behave essentially like random points. The genuine signal 
seen in the left panel of the figure is 'drowned' by the large 
number of shape-noise pair counts. The same is true for 
the cross-correlation signal depicted by the purple crosses in 
both panels. The noise-free signal is quite high, similar to 
the auto-correlation of peaks with S/M 3.4, and the error 
bars are significantly smaller than for the auto-correlation 
of the peaks with S/M 4.75 - high-S/M peaks are very 
scarce in the noise-free maps, hence the errors on their cor- 
relation function are quite large, which is the reason why we 
do not even show it in the figure. Regarding the noisy maps, 
whilst the cross-correlation function has a smaller signal-to- 
noise than the auto-correlation of the largest peaks, it can 
be used to constrain cosmology, as shown in <J4] 

3.2 The tangential-shear profiles of the peaks 

We measure the tangential-shear profiles 7t around the cen- 
tres of the peaks, as found through our filtering method. We 
use 15 angular bins, logarithmically spanning the interval 
■& e [0.3, 26] arcmin. In FigurefJ] we present the average of 
the stacked profiles, measured from 128 fields corresponding 
to the fiducial cosmology. 

The left panel shows the impact of shape noise on the 
peak profiles. The solid red pentagons and orange squares 
depict the stacked profiles of the largest peaks (S/M ^ 4.75), 
and all the peaks identified in the noise-free maps. The num- 
ber of peaks contributing to the two curves are ~ 4, 500 and 
~ 39, 000 respectively. The noise-free peaks are generated 
by lensing of large-scale structures (LSS); therefore, they are 
also present in the noisy maps, albeit with slightly different 
centre coordinates and amplitude. We match the coordinates 
of the noise- free peaks to the coordinates of the peaks in the 
noisy maps. This is done on a field-by-field basis, allowing for 
a displacement of 5 pixels around the centre of the noise- free 
peaks. The green empty squares and purple stars indicate 
the stacked profiles of the same peaks depicted by the solid 
red pentagons and orange squares, when measured from the 
noisy maps. The shape noise renders the measured profiles 
shallower on small scales, $ ^ 2 arcmin, due to the shift in 
the centre coordinates. On scales & > 2 arcmin there is no 
significant difference between the noise-free and noisy pro- 
files. The right panel shows measurements from the noisy 
maps. The <S/A/"-bins are similar to those used for the cor- 
relation function: 3.4 > S/M ^ 2.6 - green solid triangles; 
4.75 > S/M > 3.4 - blue solid squares; and S/M > 4.75 - 
red solid circles. For illustration purposes, we also include 
the stacked profiles of all peaks binned together, represented 
by the purple empty circles. The more massive the peaks, 
the higher their profiles: there is about an order of magni- 
tude difference between the average tangential shear corre- 
sponding to the highest and lowest S/M-hins,. Similar to the 
correlation function, the stacked profiles of the noisy peaks 




10 

S/N 



Figure 3. The peak function measured from the maps corre- 
sponding to the fiducial cosmological model. The results show the 
average of 128 fields of 144 deg 2 . The red points denote measure- 
ments in the absence of shape noise; the blue squares represent 
measurements from noisy maps, with an 5/A^-threshold of 2.6. 



are significantly lower than the noiseless ones. For example, 
compare the average profile of all peaks from the two pan- 
els, i.e. the orange solid squares on the left with the empty 
purple circles on the right: there is a factor of ~ 5 difference 
on small scales, and less on larger scales. This is a conse- 
quence of several factors: i) most of the \ow-S/M peaks are 
caused by shape noise and they represent random points in 
the shear maps, thereby lowering the profile average of the 
real peaks; ii) some small genuine LSS peaks are boosted by 
shape noise above the detection threshold, and so they are 
included in the right panel, but not in the left one; iii) as 
shown in the left panel, shape noise does tend to lower the 
average profiles due to the shifting of the peak centre. 



3.3 The peak function 

Figure|3] illustrates the peak function $, i.e. the number of 
peaks per unit S/M, in a field with area 144 deg 2 . We con- 
sidered 15 logarithmic S/M bins, in the interval [2.6, 14]. 
The red solid circles denote the noiseless peak function mea- 
sured from 128 fields corresponding to the fiducial cosmol- 
ogy, while the blue solid squares are the same measurement 
in the presence of shape noise. The two peak functions are 
very similar at the high-S/M end, and differ by a factor 
of at least 20 at the \ow-S/M end, where peaks generated 
by shape-noise completely dominate the genuine LSS peaks. 
Modelling the effect of shape noise on the peak function is 
the subject of a work in progress, so we shall not dwell on it 
any longer here. We merely mention that the dominance of 
shape-noise peaks at the low-S/M end of the peak function 
is in agreement with the behaviour of the peak-peak corre- 
lation function, and the tangential-shear profiles displayed 
in FiguresQ] and [2] 
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4 FISHER MATRIX CALCULATIONS 

In this section we employ the Fisher-matrix formalism to 
derive cosmological constraints from the three peak probes 
presented in the preceding section. We shall make use of the 
cosmology dependence of our simulated maps, and measure 
the peak statistics corresponding to the models varying cos- 
mological parameters around the fiducial values. These mea- 
surements can then be used to obtain Fisher predictions. 

Throughout this work we shall assume the likelihood 
function of the peak probes to be a multivariate Gaussian. 
Consider a vector of measurements of the peak-peak cor- 
relation function, peak abundance, and peak profiles m = 



{wi, . 



,<&i 



>7tj 



,7y B }, where n%, n%, n£ 



are the number of bins for the respective measurements. If 
we similarly define a vector rh of the mean of the measure- 
ments, we write the likelihood as 

1 



£(m\rh(p),C(p)) = 



(2tt) 



x exp 



W 2 [C|V2 

-(m — m) t C _1 (m — rh) 



(6) 



where the variable p indicates the dependence on the cosmo- 
logical model of the mean and covariance matrix of the mea- 
surements, in this case specified by: {n, as, f2m, w}. We have 
also introduced the total number of bins na = nj +n§ +rig T , 
as well as the total covariance matrix of the measurements 
in bins i and j: 

dj = {(rm — fhi) (rrij — ?%)>, i,j = l,n B 



(7) 



The Fisher matrix is defined as the ensemble-averaged Hes- 
sian of the logarithm of the likelihood function, 



d 2 ln£ 



(8) 



Note that we shall also use Eqs (0 and (JSj) for each probe 
individually, case in which the vectors m, fh and matrix C 
will implicitly contain only the measurements of the respec- 
tive probe, and «b will be the number of bins for that same 
probe. We shall make no further specification on this subject 
and rely on the context for clarity. 

The general exp ression of Eq. Q can be rewritten as 
jTegmark et al.ll 19971 ) 



■Fa/3 = -Tr [C 1 C, Pa C 1 C, Pfj 



m „ C m 



Pf>- 



(9) 



Here we have denoted by x tPa = dx/dp a - We shall ignore 
the first term in Eq. © for two reasons: i) it is considered to 
contribute little to the Fisher in formation - for a discussion 
on this see iTeemark et~aD l|l997h : ii) an accurate determina- 
tion of the derivatives of the covariance with respect to the 
cosmological parameters would require more realizations of 
the variational cosmologies than we currently have. There- 
fore we shall follow the standard approach to Fisher-matrix 
forecasting in the literature and ignore the trace-term in the 
above equation. 

From the Fisher matrix, one may obtain an estimate of 
the marginalized errors and covariances of the parameters: 

°"Lp<3 = [-T 7-1 ]^, (10) 
as well as the unmarginalized errors: 



-1/2 



(11) 



We now turn our attention to the estimation of each ele- 
ment contributing to the simplified expression of the Fisher 
matrix, discussed above. For each cosmological model, the 
mean auto- and cross-correlation functions are evaluated as 
an average of the correlation functions for each field, 



1 - 

>(*) = - £ «&„(*), 



(12) 



where w& is the correlation function measured from the k th 
field, and n is the number of fields considered. Similarly, we 
compute the mean of the peak function and the tangential- 
shear profiles. An unbiased, maximum-likelihood estimator 
for the covariance is 



1 - 



(rrij 



(13) 



We use this estimator to compute the fiducial-model- 
covariance matrix for each of the three probes, as well as 
for any combination of them. The estimat e of the inverse 
covariance is correc ted in the following way l|Andersonlll958l : 
lHartlap et ai1l2007n : 



c- 1 = 



JIB 



-(cy 



(14) 



n- 1 

In order to obtain low-noise estimates of the derivatives with 
respect to the cosmological parameters, we take advantage 
of the matched initial conditions of the simulations, and use 
the double-sided derivative estimator 

dfhi _ 1 ^ m'ii'Pct + Aa) — mf(p a 
~dp^ ~ n 2A^ 



Aa) 



(15) 



where Aa represents the + step in the cosmological param- 
eter p Q , e.g. Table [T] Note that whilst this approach reduces 
the cosmic variance in the derivatives, the estimated Fisher 
matrix is still noisy due to the inverse covariance estimator. 

Finally, we mention a technical point regarding the 
joined constraints presented in the next section. Since we 
evaluate the combined covariance for the three probes, we 
are concerned that its elements might be very different, 
hence introducing numerical instabilities in its inversion, and 
possibly leading to an inaccurate estimate of the inverse. To 
prevent this, we a pply the same strategy recently used in 
|Smith et alj|2012l ): instead of using the covariance matrix, 
we use the correlation matrix, defined by 



Cij j (Ji (7 j 



(16) 



where Oi = \/Cu is the rms variance. The inverse correlation 
and covariance are related by: 



We can rewrite the Fisher matrix from Eq. ([9| as: 



rm,. 



(17) 



(18) 



It is this last equation that we shall be using for our fore- 
cast, scaling the measured derivatives of the probes by the 
rms variance of the fiducial functions, and employing the 
inverse correlation matrix, instead of the inverse covari- 
ance. Figure[4] depicts the derivatives of the auto- and cross- 
correlation function with respect to the cosmological param- 
eters that we investigate in this work - left and right panels, 
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10 100 10 100 

i3[arcmin] i9[arcmiri] 

Figure 4. The derivatives of the peak correlation function measured from the noisy maps. The left panel shows the auto-correlation 
function ujp p for peaks with S/M 5= 4.75, while in the right one we present the cross-correlation u) pp between two bins with low and high 
S/M. 




1 10 1 10 

i9 [arcmin] T?[arcmin] 



Figure 5. The derivatives of the tangential-shear profiles -jt measured from the noisy maps. The left panel corresponds to the peaks 
detected with S/M > 2.6, while in the right one, peaks have S/M > 4.75. 



respectively. The auto-correlation derivatives are quite noisy 
- primarily due to the relatively small number of peaks with 
S/M ^ 4.75, a characteristic also visible in Figure[T] the 
noise being largest on small scales, $ < 20 arcmin. The rim- 
derivative has the largest amplitude, while w has the small- 
est. All derivatives tend to on larger scales, signifying that 
there is no information in the peak-peak correlation function 
for $ > 2 deg. Compared to the auto-correlation function, 
the derivatives of the cross-correlation are less noisy, with 
the exception of n s . They are also significantly smaller, just 
like the cross-correlation signal is much smaller than the 
auto-correlation one. In both panels, the fi m -derivative is 



negative for i? ^ 10 arcmin, whilst for the other parameters, 
the derivatives cross the 0-line a few times. 

Figured] presents the shear-profile derivatives with re- 
spect to the cosmological parameters: the left panel corre- 
sponds to peaks with S/M 2.6, while the right one is for 
S/M 4.75, hence the larger measurement noise present 
there. In both panels, the derivatives decrease to for in- 
creasing angular scales, e.g. $ ~ 20 arcmin. In the left panel, 
all derivatives are positive, with the only exception of w, 
which has a barely-noticeable transition from slightly neg- 
ative to slightly positive. The fl m -derivative is the largest, 
while the w-one is the smallest. For the higher-<S/A/" peaks, 
the derivatives are less smooth, but the general trends are 
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Figure 6. The derivatives of the peak function <E> measured from 
the noisy maps and scaled by the peak function corresponding to 
the fiducial cosmological model. 



preserved: f2 m and as still have the highest derivatives, fol- 
lowed by n s and w. They all seem to peak on the scales 
0.8—1 arcmin, which suggests that a significant amount of 
information arises from these scales. We shall mostly be us- 
ing the larger peaks for our constraints. 

Figure[6] depicts the derivatives of the peak function 
with respect to the cosmological parameters, scaled by the 
peak abundance corresponding to the fiducial cosmology, as 
a function of S/N. Again fi m has the largest derivative, fol- 
lowed closely by erg - both of them positive. The derivatives 
with respect to n B and w are much smaller and less smooth. 

Since in this study we use numerical derivatives to es- 
timate the Fisher errors, as opposed to analytical ones, we 
need to ask what the resulting uncertainty on the Fisher 
errors is. This is particularly needed for our noisiest mea- 
surements, i.e. for the correlation function. We address this 
question in the appendix SjX] and refer the interested reader 
to that discussion. 

One point of interest - useful mainly to develop one's 
intuition - is to explore which are the highest contribut- 
ing scales to the Fisher information. Naturally, we expect a 
combination of scales to be the most effective at constraining 
the cosmological model; the most common approach to de- 
termine such a combination is to diagonalize the covariance 
matrix and determine its eigenvalues: the eigenvectors cor- 
responding to the largest eigenvalues will represent the most 
constraining linear combinations of scales. To see if there is 
a particular scale which helps to significantly reduce the er- 
rors, we shall take two approaches. First, we compute the 
cumulative Fisher matrix, i.e. we monitor how the errors 
evolve when bins are included one-by-one in the calculation 
of Eq. ©. We choose as starting bins those containing mea- 
surements expected to be easier to perform in real surveys. 
For example, for the peak counts, we start with the highest- 
S/N bins, and then add progressively sm&Wev-S/N bins. For 
the profiles and correlation functions, we start with small- 
scale measurements, and progressively include large angular 



scales. The results of this exercise are presented in the left 
panels of Figure]?] from top to bottom, we show the cumula- 
tive unmarginalized Fisher errors arising from cjp p , yr, and 
$, scaled by the final errors - when all bins are included 
- as a function of bins. For the auto-correlation function, 
the information seems to saturate at ~ 40 arcmin: the in- 
clusion of larger-scale measurements does not further im- 
prove the constraints. For the profiles, the saturation occurs 
at ~ 10 arcmin, while for the peak counts, the inclusion of 
peaks with S/N < 4 still improves the constraint on w, 
though not o n the other paramete rs. This is in line with our 
findings from iMarian et all (|2012t >. 

The second approach involves two approximations: i) 
for each probe we neglect the covariance of the parameters, 
i.e. we consider only unmarginalized errors - this was also 
done for the first approach, ii) we disregard the off-diagonal 
elements of the covariance matrix of each probe. The Fisher 
matrix can then be written as 



(19) 



We can quantify the contribution of each bin j to the total 
Fisher constraints through the ratio: 



T 3 



ICj; 



jCii 



(20) 



By construction, this ratio works best as a contributing- 
scales indicator for those probes which have a covariance 
matrix with as few off-diagonal elements as possible. In 
our case, this is the auto-correlation function, as can be 
seen from Figure|8] The rms of the above ratio is shown 
in the right panels of Figure[7J ordered in the same way as 
the left ones. For the correlation function, the most con- 
tributing scales to the as- and f2 m -constraints are ~ 30 and 
~ 40 arcmin respectively. This is in agreement with the left 
panel of the figure: the errors on those two parameters are 
lowered significantly when the respective bins are included 
in the Fisher estimates. The constraints from stacked pro- 
files receive contributions from all scales below 10 arcmin, 
while for the counts the smaller bins (S/N < 5) tighten the 
errors - though notice that the agreement between the left 
and right panel is not so good as for the correlation func- 
tion, due to the high correlations between the low-S/N bins. 
These correlations are included in the left panel, but not in 
the right one. 



5 COSMOLOGICAL CONSTRAINTS 

Using all the ingredients described in section iJU we now pro- 
ceed to constraining the cosmological model using the three 
peak probes individually, as well as in combination. To this 
avail, we shall employ only noisy maps, and consider the 
auto-correlation wjj p of peaks with S/N ^ 4.75, the cross- 
correlation ojpp of peaks with S/N > 4.75 and 3.4 J: S/N ^ 
2.6, the stacked profiles 7t of peaks with S/N ^ 4.75, and 
the abundance of peaks $ with S/N > 2.6. The Fisher ma- 
trix errors are estimated using the covariance on the mean 
of the fiducial model for the three probes of interest to us 
here. Thus, our estimated covariance from Eq. (|13p . which 
corresponds to an area of 144 deg 2 , is rescaled to correspond 
to an area 128 times larger, i.e. ~ 18000 deg 2 . Together with 
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Figure 7. Left panels: Cumulative unmarginalized Fisher errors for <^p P , 7t> an< i ^ _ top, middle, and bottom panels respectively. Bins 
are added progressively, from small to large scales for the angular bins, and from large to small S/J\f. The cumulative errors are scaled 
by the final error, i.e. the error obtained when all bins have been included. For the auto-correlation and the profiles, only peaks with 
S/Af ^ 4.75 are included. Right panels: The contribution of different scales to the Fisher information, estimated as the rms of Eq. 1 1201 1. 
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Figure 8. Total cross-correlation matrix of the peak abundance 4>, the peak-peak correlation function ^p P , and the stacked peak profiles 
fr- w p P an d 7T are evaluated for peaks with S/N > 4.75, while $ contains the peaks with S/N > 2.6. From left-right and down-up, 
we plot the correlations of <I>, ^p P , 7t, with the S/Af and angular scales (measured in arcmin) increasing in the same directions. The 
correlation matrix is computed from 128 fields of 12 x 12deg 2 corresponding to the fiducial cosmology and is rescaled to match a sky 
coverage of ~ 18000 deg 2 . 



the survey specifications given in section [JU this makes our 
study representative for two future surveys, LSST and Euclid. 

Figured] presents the cross-correlation matrix r of these 
probes. By far, 7t has the strongest correlation coefficient 
of the three: ~ 0.7 on scales 2 — 20 arcmin. For the peak 
function, the lavf-S/Af bins are the most correlated ~ 0.5 
for S /Af < 5. This was already established in our earlier 
work (Ma rian et alj[2012h . and it can be explained through 
the better-known behaviour of halos: small-mass halos are 
sample- variance dominated, while t he large and rare ha- 
los follow the Poisson distribution (|Hu fc Kravtsovl 120031 ; 
ISmith fc Marianll20lih . Note however how the smallest-S/A/" 
bins in Figure|8] seem to be completely uncorrelated: this is 
most likely due to the overwhelming number of shape-noise 
peaks, which are random, unclustered, and therefore uncor- 
related. o»p P displays the smallest correlation coefficient of 
the three, ~ 0.3 — 0.4 on the scales 20 — 60 arcmin, with 
weaker correlations on smaller scales. We further note the 
weak cross-correlation of ojp P and as well as ojp P and 7t- 
There is a visible cross-correlation of $ and yr, of ~ 0.3 for 
peaks with S/N > 7. This is most likely due to the stacked 
profiles being dominated by the most massive peaks, which 
also dominate the high-S/N end of the peak function. Ta- 



ble [2] presents the unmarginalized and marginalized 1-a er- 
rors resulting from the three peak probes. Each probe taken 
by itself, the abundance of peaks has the greatest constrain- 
ing power, followed by the profiles, and then by the corre- 
lation function. Regarding the latter, we note that ojp p and 
cdpp yield very similar constraints, the auto-correlation be- 
ing more effective for w and £7 m - a reduction by factors 
of ~ 2 and ~ 1.5 respectively in these errors, compared 
to the cross-correlation. However, when combined with the 
other two probes, there is little difference between Ljp P and 
ojpp. The greatest benefit to adding the correlation func- 
tion or the profiles to the abundance of peaks concerns the 
time-independent equation-of-state for dark energy: after 
marginalizing over the other parameters, the errors on w re- 
sulting from 3> and Up p taken individually are similar, while 
the profiles seem to yield a constraint tighter by a factor of 
~ 1.7. When all three probes are combined, the constraints 
on Q m , as, n s improve by a factor of ~ 1.5 — 2 compared to 
using $ alone, while for w the improvement is ~ 2.5. Lastly, 
combining $ and 7t is almost as efficient as using all three 
probes: the contribution of the correlation functions to re- 
ducing the errors is negligible, if both the abundance and 
the profiles are used. 
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Figure 9. 2-<r Fisher ellipses for the peak probes. The auto-correlation function is depicted by the green dot-dashed ellipse; the red 
dashed ellipse represents the profiles, while the dotted blue line corresponds to the peak abundance. The combination of the three is 
given by the black ellipse. 



Table 2. Fisher-matrix constraints using Wp p , u)^ 
parameters, f2 m = 0.25, erg = 0.8, w = —1, n B = 1. 



7t, and <E>. All errors are presented as percentages of the fiducial values of the 
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This is further detailed by Figure[9] which shows the 
forecasted likelihood contours, at 95% confidence level. The 
largest ellipses - green, dot-dashed - represent the ujp P con- 
straints; for many of the parameter combinations, the red, 
dashed ellipses depicting jt are almost perpendicularly ori- 
ented relative to the ojp P ellipses, and also quite tilted with 
respect to the <3? ellipses, shown by the blue dotted lines. 
The greatest gain in combining all three probes is for the 



case of w, where the resulting ellipse is smaller by a factor 
of ~ 3 compared to the $ ellipse. 

We compare our Fisher constraints to those of 
Iffilbert et ail (|2012l ). derived for an identical survey sce- 
nario from the same simulated WL maps. That analysis com- 
bines the shear correlation functions, the third moment of 
the aperture mass, and the peak abundance, and accounts 
also for local primordial non-Gaussianity. The percentage 
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Figure 10. 2-cr Fisher ellipses for the peak probes, with a PLANOK-like prior. The colors and symbols are the same as in the previous 
figure. The prior has been added to each peak probe separately. 



unmarginalized combined 1-a errors found in iHilbert et al.l 
(|2012T) for {n m , a 8 , w, n s } were {0.2, 0.15, 0.9, 0.4}. The 
marginalized errors were {0.7, 0.4, 2.6, 0.7}. Their con- 
straints are tighter, particularly for fi m , but the combination 
of peak probes performs nonetheless very well. 

Figure[lO] presents the 95%-likelihood contours of the 
three probes, after we have added cosmic microwave back- 
ground (CMB) information resulting from an experiment 
similar to PLANCK. For the latter, we assume that the CMB 
temperature and polarization spectra can constrain 9 pa- 
rameters: the dark energy equation-of-state parameters wo 
and w a ', the density parameter for dark energy Ode; the 
CDM and baryon density parameters scaled by the square 
of the dimensionless Hubble parameter cjcdm = ^cdm h 2 
and uih = flbh 2 (h = _ffo/[100 km s _1 Mpc -1 ]); the primor- 
dial spectral index of scalar perturbations n s ; the primor- 
dial amplitude of scalar perturbations A 3 ; the running of 
the spectral index a; and the optical depth to the last scat- 
tering surfac e r. We comput e the CMB Fisher matrix as 
described bv lEisenstein et all (1 19991 ): 

= ^Cov- 1 [C t)X ,Ci, Y ] ^ , (21) 



l X,Y 



dp a 



where {X, Y} e {TT, EE, TE, BB}, where C ; ,tt is the tem- 
perature power spectrum, C^ee is the E-mode polariza- 
tion power spectrum, C;,te is the temperature-E-mode po- 
larization cross-power spectrum, and C;,bb is the B-mode 



polarization power spectrum. The assumed sky coverage 
is /sky = 0.8 In order to make the CMB Fisher matrix 
compatible with our parameters, we rotate it to a new set 
q T = {w ,w a ,Q. m ,h, f b ,r,n s ,a$,a} , where for us w = w. 
We marginalize over the 5 parameters absent from our 
analysis, and then add the resulting 4x4 matrix to the 
Fisher matrix for each of the probes, and their combina- 
tion. This way of implementing the PLANCK prior follows 
the same steps taken in o ur previous studies (|Marian et al.l 
120121 ; IHilbert et all [20121 ) . FigurefTDI shows that combining 
the peak probes with the CMB dramatically alters all the el- 
lipses involving n s : they are shrunken and almost completely 
aligned. This is hardly surprising, as the CMB constrains n s 
much better than any WL probe, so the ellipses are then 
dominated by the CMB information. For the other param- 
eters, the changes are small, and there is still a significant 
gain to be obtained by combining the peak probes, compared 
to using only one of them, for instance the abundance. 



6 SUMMARY AND CONCLUSION 

In this paper we have investigated cosmological constraints 
from WL peaks, beyond using just counts. We have mea- 
sured for the first time the peak-peak correlation function, 
and estimated its constraining power on the following cos- 
mological parameters: {f2 m , ag, w, n s }. We have also mea- 
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sured the tangential-shear profiles of the peaks, and explored 
their cosmological utility. Whilst both the correlation func- 
tion and stacked profiles of galaxies/clusters are standard 
tools for cosmology, our study applies them to shear-selected 
objects, detected only through their gravitational lensing 
effects. We employed mostly numerical methods, perform- 
ing measurements of the peak abundance, peak profiles, and 
peak-peak correlation functions in simulated WL maps. To 
estimate the cosmological constraints, we used the Fisher- 
matrix formalism, computing the derivatives with respect 
to cosmological parameters and the combined covariance of 
the three peak probes from the simulated maps. The lat- 
ter have been generated with ray-tracing through a large 
suite of iV-body simulations, varying the cosmological model 
around some fiducial values. Given the distribution of source 
galaxies and the level of shape noise considered, our simu- 
lated maps are relevant for future large WL missions, such 
as Euclid and lsst. 

To measure the correl ation function of peaks, we used 
the lLandv fc Szalavl (|l993h estimator, known for two essen- 
tial properties: i) minimum variance; ii) lack of bias due 
to the unknown mean density of the measured objects. We 
have shown that the correlation function has its maximum 
on scales # < 10 arcmin, with the position and height of 
the latter depending on the S/Af of the considered peaks. 
For angular scales i9 ^ 8 arcmin, the correlation function of 
'large' peaks (S/Af ^ 4.75) has a greater amplitude than 
that of 'medium' and 'small' peaks (4.75 > S/Af ^ 3.4 and 
S/Af < 3.4), mirroring the clustering behaviour of dark- 
matter halos. We have shown that the low-S/Af end of the 
peak function is completely dominated by shape-noise peaks 
- only ~ 5% of peaks with S/Af ~ 3 are genuine LSS peaks. 
This is the reason why the correlation function of the small 
peaks quickly declines to 0: shape-noise peaks are random, 
and therefore uncorrelated. On the contrary, the correlation 
function of the peaks measured from noise-free maps is quite 
high, even with all peaks, small and large, binned together: 
this is just a statement that peaks due genuinely to projected 
LSS are quite clustered. The cross-correlation of small- and 
large-iS/A/" peaks measured from noisy maps is also smaller 
than its noise-free counterpart, but it can still be used for 
cosmological constraints. 

We have measured the tangential-shear profiles of the 
peaks using their WL centre, as detected by our hierarchi- 
cal method. We found that the larger the S/Af of the peaks, 
the higher their average profile. A comparison of the results 
from noise-free and noisy maps showed that: i) the ratio of 
the average profile of high-iS/TV" peaks to the average profile 
of all peaks combined together is ~ 2 for noise-free maps, 
and ~ 5 for noisy maps; ii) the profile average of all noisy 
peaks binned together is a factor of ~ 5 smaller than that of 
the noiseless profiles; iii) the profiles of genuine LSS peaks 
are shallower on scales ■& < 2 arcmin when the peaks are de- 
tected and measured from noisy maps, rather than noise-free 
maps - this owes to a coordinate shift of the peak centres in 
the presence of shape noise. The first two findings are mostly 
due to the overwhelming number of shape-noise peaks which 
dominate the \cw-S/Af end of the peak function: measure- 
ments of the tangential-shear profiles around random points 
in the map bring down the average of the stacked profiles of 
genuine LSS structures. 

The Fisher-matrix analysis that we performed revealed 



the forecasted errors on the cosmological parameters from 
the auto- and cross-correlation functions of the peaks, as 
well as from the stacked profiles to be larger by at least 
a factor of 2 than those obtained using the peak func- 
tion. Nevertheless, combining the peak function with any 
of the other three probes reduces the errors significantly, 
particularly for w. The ratio of the marginalized 1-cr con- 
straints from $ alone to the constraints from uj pp + $ is 
{1.36, 1.43, 1.5, 1.3}, corresponding to {f2 m , erg, w, n s }; for 
7t + $ the ratio is {1.6, 1.7, 2.3, 1.6}; and for tj pp + j T + 4>, 
the result is {1.7, 2, 2.5, 2}. Therefore, one of the main con- 
clusions of this study is that future WL surveys should 
use peak statistics beyond the 1-point function. Note that 
our cjp p + 7t + $ constraints are qu ite competitive with 
those derived by iHilbert et alj (|2012T ) from combining the 
shear correlation functions, the third moment of the aper- 
ture mass, and the peak abundance. For the parameters that 
we have investigated, the stacked tangential-shear profiles 
seem slightly more constraining than the correlation func- 
tion. However, if other parameters were to be included in 
the analysis, particularly those quantifying primordial non- 
Gaussianities, the peak-peak correlation function might be 
more powerful. This however, is the subject of an ongoing 
work, to be presented in the near future. 

We have briefly examined the S/Af and angular scales 
contributing most to the Fisher information. The cosmol- 
ogy constraints from the peak-peak correlation function 
benefit significantly from measurements at scales i? ~ 
20 — 40 arcmin, and saturate for ■& > 50 arcmin. For 
the tangential-shear profiles, measurements at 1? ~ 1 — 
10 arcmin are important, and the constraints saturate for 
■& > 10 arcmin. The peak abundance benefits from the in- 
clusion of smalLS/TV bins. 

In this study, we have not exhaustively investigated the 
properties of the clustering of the peaks. Measuring the cor- 
relation function was only the first step of such a study. 
It would also be interesting to examine the biasing of the 
peaks with respect to the shear/convergence field, and also 
to compare the clustering of the peaks with that of the halos 
in the simulations. Our goal for the near-future is to deter- 
mine the full constraining power of the shear-selected WL 
peaks on an expanded set o f cosmological param eters, i.e. to 
extend the recent study of IHilbert et al.l l|2012f) to include 
the 2-point statistics of the peaks, as well as tomographic 
techniques. 
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APPENDIX A: ESTIMATING THE ERROR ON 
THE FISHER CONSTRAINTS 

In Fisher forecasts as well as real-data analysis, the deriva- 
tives of observables with respect to cosmological parameters 
are usually estimated analytically. However, due to the lack 
of analytical predictions for peak observables, in this study 
we employ numerical derivatives. As discussed in §3] the 
latter are estimated as double-sided derivatives (Eq. (|15|) ). 
which reduces the impact of cosmic variance and the de- 
pendence on the step Ap a in the cosmological parameters 
around the fiducial value. One natural question is: what is 
the impact of the uncertainties in the measured derivatives 
on the resulting Fisher constraints? We shall address this is- 
sue in an approximate manner, just for the unmarginalized 
errors. 

The unmarginalized Fisher constraints are computed 
from Eqs (|ll|l and (|18[) . Their variance can be written as 

Var[a PQ ] = i J^ Q 3 Var^] . ( Al) 

Since we are interested in the impact of the uncertainties 
in the derivative estimates on the Fisher matrix, we ignore 
the uncertainties in the covariance matrix of the probes. 
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In any case, we would not be able to reliably estimate the 
latter from just 128 realizations, so this approximation is 
quite necessary. Note that the impact of uncertainties in 
the covariance mat rix on Fisher constr aints was explored in 
the recent study of lTavlor etafl (|2012T ). The variance of the 
Fisher-matrix elements can then be written as 

"b / r -i \ 2 
VarLF^] = V { m j,j>sVar[mi, Pot ] 

+ r "'i,p Q Var[mj, P/3 ] + 2 m i;Pcl m J)P( , Cov[m, !Po mj, P|J ]}. (A2) 

The rms of the variance of the derivatives represents the 
error bars shown in Figures[4l [5] and [6] The covariance in 
the above equation must be estimated separately from the 
data. Since we have only 64 realizations corresponding to 
each of the variational cosmologies, this estimate is bound 
to be very noisy. Whilst fn can be any of the three peak ob- 
servables, we actually perform the above calculation for the 
correlation function alone. The derivatives of the latter have 
the largest error bars and they also float around the line, 
which could yield spurious Fisher information. Here we men- 
tion that increasing the number of bins to 12 and 15 did not 
significantly change our parameter constraints for u>p p . The 
uncertainties in the unmarginalized l-cr Fisher constraints 
for Wp P , expressed as percentages of the values given in Ta- 
bleware for {ft m , as, w, n s ): {26%, 31%, 17%, 26%} respec- 
tively - if we ignore the noisy covariance term in Eq. (|A2[) ; 
and {38%, 43%, 24%, 36%} - if we do take into account the 
covariance term. 

We leave to a future work a more complete treatment 
of the Fisher-matrix constraints from numerical measure- 
ments. We would like to analytically model the peak probes 
so as not to be vulnerable to numerical effects when com- 
puting their derivatives with respect to cosmological param- 
eters. 



© 0000 RAS, MNRAS 000, 000-000 



